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ABSTRACT 



Two distinct scenarios for the origin of the ~ 4 x 10^ Mq of dust observed 
in the high-redshift {z = 6.4) quasar J1148+5251 have been proposed. The first 
assumes that this galaxy is much younger than the age of the universe at that 
epoch so that only supernovae (SNe) could have produced this dust. The second 
scenario assumes a significantly older galactic age, so that the dust could have 
formed in lower-mass asymptotic giant branch (AGB) stars. Presenting new in- 
tegral solutions for the chemical evolution of metals and dust in galaxies, we offer 
a critical evaluation of these two scenarios, and observational consequences that 
can discriminate between the two. We show that AGB stars can produce the 
inferred mass of dust in this object, however, the final mass of surviving dust de- 
pends on the galaxy's star formation history (SFH). In general supernovae cannot 
produce the observed amount of dust unless the average SN event creates over 
~ 2 Mq of dust in its ejecta. However, special SFHs can be constructed in which 
SNe can produce the inferred dust mass with a reasonable average dust yield of 
~ 0.15 Mq. The two scenarios propose different origins for the galaxy's spectral 
energy distribution, different star formation efficiencies and stellar masses, and 
consequently different comoving number densities of J1148+5251-type hyperlu- 
minous infrared (IR) objects. The detection of diagnostic mid-IR fine structure 
lines, and more complete surveys determining the comoving number density of 
these objects can discriminate between the two scenarios. 

Subject headings: galaxies: evolution, high-redshift, starburst - quasars: individ- 
ual: (SDSS Jl 14816.64+525150.3) - infrared: galaxies 
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INTRODUCTION 



Determining the origin of the massive amount of dust present in SDSS J1148+5251 
(hereafter J1148+5251), a hyper-luminous quasar at z = 6.4 presents a special challenge. 
The dust mass inferred from far-infrared and submillimeter observations is ~ (1 — 5) x 10^ Mq 
depending on dust composition (Dwek et al. 20 07, and references th erein), about ten times 
larger than the mass of dust in the Milky Way ( ISodroski et al.lll997l ). The total far-infrared 
(IR) luminosity is about 2 x 10^'^ Lq, giving a star formation rate (SFR) of ~ 3400 M0 yr~^ 
for a Salpeter initial mass function (IMF) (IKennicuttl Il998al ) . Since the universe was only 
~ 890 Myr old at that redshift, and the galaxy perhaps much younger, it has been suggested 
that only c ore-collapse supernovae (CCSNe) can produce the observed arnount of dust in 



2004; Rho et al 



20081 



Sibthorpe et al 



this object (Morgan fc Edmunds 2003 : Dunne et al.ll2003l : iNozawa et al.ll2003l : iMaiolino et al. 



2OO9I ). Another argument made in favor of SNe as 



the most important dust sources in the early universe, was based on observations of unmixed 
ejecta of young remnants, primarily Cas A. The inferred dust mass in this remnant, ranging 
from ~ 0.1 — 0.15 Mq, was considered sufficiently large for SNe to account for the obs erved 
dust mass in J1148+5251 jRho et aLllioOsI : Isibthorpe et~al1 12OO9I: ISarlow et allboioh . We 
consider the 



1 Mq of dust mass claimed by 



Dunne et all (2003 



20091 



to have formed in 

Cas A as an unreasonable amount of dust. The tot al ejecta mass i s 2 — 4 Mq, of which 



only ~ 0.2 Mq consists of condensible elements [e.g. iNozawa et al.l ( 120101 ). and references 
therein] . 

None of these assertions, that SNe are significant dust sources in the early universe, 
were substantiated by detailed calculations. Specifically, the role of supernovae as destroyers 
of dust during the remnant phase of their evolution was completely ignored. Furthermore, 
these claims tacitly assume that the epoch of intense star formation spanned the entire age 
of the galaxy, which requires an excessively large reservoir of interstellar gas. Finally, recent 
calculations show that on ly ~ 0.15 Mq of dust is created in a Population III 20 Mq SN 
Jcherchneff fc Dwekl lioiol ) . which may be a typical yield for Pop II SNe of similar masses. 
All these considerations point to the need of more detailed calculations to ascertai n the role of 



SNe as dust sources and sinks in the early universe. In such recent calculations, (IDwek et al. 



20071 ) explored the combined effects of dust formation and destruction on the net amount of 
dust produced in J1148-I-5251. They found that SNe cannot be important sources of dust in 
the early universe unless they produce significantly more, and destroy significantly less dust 
than implied by current observations or theoretical calculations. 



In light of these difficulties, IValiante et al.l (120091 ) suggested that AGB stars could be the 
source of dust in J1148-I-5251. Their model was motivated by the numerical simulations of 
the formation and growth of this galaxy through a series of successive mergers that resulted 
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in repeated intense bursts of star formation [ (ILi et al.l 120071 ). hereafter Li07]. Associated 
with these mergers is the growth of its central black hole (BH) which at z ^ 6.4 has reached 
a mass of about ~ 10^ Mq. In this scenario, star formation commenced at z ~ 15, when the 
universe was merely 250 Myr old. Consequently, the progenitors of the more numerous and 
efficient dust producing AGB stars had time to evolve off the main sequence (see Table 1), 
and produce the observed mass of dust in J1148+5251. 

This paper takes a critical look at these two proposed dust formation scenarios, hereafter 
referred to as the SN and AGB scenarios. We start by introducing some basic definitions 
of the various quantities that govern the chemical evolution of galaxies, and present new 
integral solutions for the chemical evolution of their elemental and dust content (Section 2). 
In Section 3 we discuss the yields of the main stellar dust sources. These include the explosive 
SN ejecta, the fast winds created during the Wolf-Rayet (WR) stage of the evolution of 
stars with masses in excess of ~ 40 Mq, and the quiescent winds from AGB stars. In 
Section 4 we apply the model to follow the chemical evolution of Jl 148+5251. We calculate 
the contribution of SNe and AGB stars to the production of dust in this galaxy, and the 
evolution of its stellar mass and luminosity. The simulated SFH of Li07 leading to the 
formation of J1148+5251 is not unique, and in Section 5 we examine alternate scenarios for 
the formation and survival of dust by AGB stars and SNe. Our results show that both, 
the AGB and the SN scenarios, are still viable for producing the mass of dust observed in 
J1148+5251, each scenario with its own problems and limitations. We therefore present in 
Section 6 observational tests, including studies of the spectral energy distribution (SED), 
the rarity of J1148+5251-type objects, and the inferred SFR, that can discriminate between 
the two. In Section 7 we explore two additional, non-stellar, dust sources: molecular clouds, 
which have always been considered as an environment for the growth and processing of 



inter s tellar dust grains in the Galaxv (ISnow 



19891 : lDweklll998l : iGreenberg fc Lilll999l : IZhukovska et al. 



1975 



Dwe 



c fc Scald Il980t iLiffman fc Clayton 



20081 ) , and AGN winds, which were 



recently suggested a s potential producers of interstellar dust in quasars (lElvis et al.l 12002 



Maiolino et al.ll2006[ ). The results of the paper are briefiy summarized in Section 8. 



Throughout this paper we adopt a fiat ACDM cosmology, with a baryonic density 
parameter Qi, = 0.044, a total matter (dark+baryonic) density parameter of Qm = 0.27, 
a vacuum energy d ensity = 0.73, and a Hubble constant of Hq = 70 km s~^ Mpc~^ 
jSpergel et al.ll2007l ). 



-4- 



2. EQUATIONS FOR THE EVOLUTION OF DUST 

2.1. Basic Definitions 

We define the stellar IMF, 0(m), so that 0(m)(im, is the number of stars with masses 
between m and m + dm, and normalize it to unity in the {mi, m„} mass interval, where mi, 
and mu are, respectively, the lower and upper mass limits of the IMF. The IMF-averaged 
stellar mass, (m), is then: 

{m) = m (f){m) dm , (1) 

J nil 

The SFR, ip{t), is the mass of stars formed per unit time, and is related to the stellar 
birthrate, B{t), by: 

Bit) = (2) 

(m) 

We assume that all stars in the [mi, mj\ mass range end their life quiescently, whereas all 
stars with masses m^ < m < m^ become CCSNe, where = 8 Mq, is their lower mass 
cut. 

A useful quantity is m^, the mass of all stars born per SN event, given by: 

rnriu 

m^ = (m) / / (f){m) dm (3) 

The SN rate, Rsn, is then given by: 

i?SN(t) = B{t) / 0(m) dm = ^ , (4) 

In all our calculations we will use a mass-heavy ISM characterized by a power law: 
0(m) ~ m~" in the {mi, mu} = {1 Mq, 100 Mq} mass interval, with a = 2.35. For this 
IMF we get: 

(m) = 3.1 M0 and m^ = 53.0 M© (5) 
The choice of this mass-heavy IMF is motivated by studies s uggesting that star forin ation 



in the early universe was biased towards more massive stars f lBromm Sz Larson! |2004| ). For 



comparison, a Salpeter IMF gives values of (m) ^ 0.35 Mq, and ^ 135 Mq. 



2.2. Equations for the Evolution of Dust or Elements 



We first describe the equations for the evolution of the mass of dust in the ISM. Let 
Mg{t) be the mass of the gas in the ISM at a given time t. The evolution of M^it), the mass 
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of dust in the ISM, is governed by the equation: 

at Td[t) 
where Zd{t) is the dust-to-gas mass ratio: 



dMdjt) 
dt 



(6) 



in f /out 



Zdit) 



(7) 



The first term in Equation ([6]) represents the rate at which the dust is removed from the ISM 
by star formation; the second is the rate at which dust is destroyed by sputtering or grain- 
grain colhsions in SN blast waves. The parameter Td{t) is the timescale for the combined 
effect of these processes; the third term, S{t), is a source function representing the rate of 
dust formation in the different astrophysical environments; and the fourth term represents 
the rate of increase/decrease in the dust mass as a result of infall/outflow from the galaxy. 



The timescale for grain destruction in supernova shock waves is given by (IDwek &: Scalo 



1980): 



Td = 



{nid) i?sN 

The parameter (m^) is the total mass of the refractory elements, initially locked up in dust, 
which are returned back to the gas phase of the ISM t h rough out the evolution of a single 
supernova remnant (SNR). For the MW galaxy, iJonesI (120041 ) found dust lifetimes of 400 
and 600 Myr for silicate and carbon grains, respectively. We adopt here an average lifetime 
of Td ~ 500 Myr. This average should be time dependent since the dust composition, 
specifically the relative abundance of silicate and carbon dust, evolves with time. However, 
such distinction is not warranted in the present investigation because of our current lack of 
knowledge of the nature of the dust in high-redshift gal axies. Adopting a total dust mass 
Md^3x W Mq, and a Galactic SN rate i?sN ~ 0.02 jPiehl et al.ll2006l ). gives a value of 
about 3 Mq for the grain destruction efficiency in the Galaxy. Since (m^) scales linearly 
with the mass of ISM dust, we can normalize it to its Galactic value: 



{rrid) 



Zd{t) 



-'d,mw 



{md)r 



where (rrid)^^ = 3 Mq is the Milky Way value of (m^), an d Zd^mw 
dust-to-gas mass ratio in the Milky Way (jZubko et al.l 120041 ). 
The grain destruction rate can then be written as: 



(9) 



0.007 is the average 



Mdjt) 
Tdit) 



Zdit) 



-RsN 



■'d,mw 



(10) 
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The source function S{t) represents the rate at which the dust mass in the ISM increases 
by nucleation in the different dust sources, or by growth in the ISM, which is an important 
process responsible for the differential depletions in its different phases. Considering only 
stellar sources, the source function can be written as: 



Sit) 



m{t) 



+ 

+ 



''wr 



(m) 
(m) 
(m) 



A 



la 



f 




JO 


(m) 



Yd^agnifn) (f){rn) dm (AGB stars) 

yd,wr{^) (pijn) dm (WR stars) 

Yd^snim) (j){m) dm (SNe) 

Yd,Iafla{T)dT (SNIa) 



'ID 



where the different terms represent the net contribution of the different stellar sources to 
the dust mass in the ISM. For each stellar mass, the SFR, ^/'(t'), is calculated at the epoch 
t' = t — Tms{m), where Tmsi"^) is the main sequence lifetime of a star of mass m (Table 1). For 
AGB stars, the lower limit of the integral, m{t), is given by: maxjm;, m'{t)}, where m'{t) is 
the mass of a star with a main sequence lifetime t. The lower limit, m^^r is the limiting mass 
above which stars undergo extensive mass loss and become WR stars. Yd^agb{m), Ya^wrifn), 
Yd^sniiTT^) , and l^./a are the dust yields in the quiescent winds of AGB and WR stars, and 
in the explosive ejecta of core-collapse and Type la SNe, respectively. In the last term, Aia 
is a normalization constant that can be determined from the ratio between the observed 
frequency of Type II and Type la SNe, and fiai^) is the distributio n function of the de lay 
time, r, between the birth of the stellar system and the SNIa event fjGreggidl2005l . |2010| ). 



Finally, the yields in eq. flTTl) . except the SNIa yields, depend on the initial stellar metal- 
licity. This implicit time dependence is suppressed in the equations, and its implementation 
in the solution is described in Section 2.4 below. 

For stellar masses above ~ 8 Mq, Tms{f^) « t for most times of interest, so that the 
yield from SNe and WR stars can be set equal to RsNYd^sn, and RwRYd,wr, respectively. 
RwR is the rate of WR stars, which can be defined similarly to the SN rate [see eq. (j4])] with 

replaced by m^r- The IMF- averaged yields of SNe and WR stars are given by: 



Yd,x = I Yd^x{m) (f){m) dm/ / (f){m) dm 

' mx J mx 



(12) 



where X = sn or wr. 
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The term describing the rate of change in the dust mass caused by the exchange of gas 
with the intergalactic medium depends on whether the gas is flowing in or out of the galaxy. 
In an outflow, the dust-to-gas mass ratio in the outflowing gas is about equal to the galaxy's 
value, Zd{t), and the term is given by: 



dM,{t) 
dt 



Zd{t) 



out 



dMg{t) 

Jt 



for outflows 



(13) 



out 



In the case of infall, the dust-to-gas mass ratio of the infalling gas, Z^J^-^ (t), is an independent 
variable and the term becomes: 
'dMdit) 



dt 



J inf 



dMg{t) 
Jt 



for infalls 



(14) 



inf 



2.3. The Evolution of the Gas 



The evolution of the ISM gas is given by the equation: 



dMg{t) 
dt 



-m + 



nil 



(m)) 



-{l-Rej)^P{t)± 



(m) 

dMg{t) 
Jt 



Mej{m) (f){m) dm ± 



dMg{t) 
dt 



(15) 



inf /out 



- inf /out 

where the first term in the top line represents the conversion rate of the ISM mass into stars, 
the second term represents the rate of stellar mass loss, where M, 



m 



is the total mass 

returned back to the ISM by a star of initial mass m, and the third term represents the rate 
of change in the ISM mass as a result of gas infall or outflow. The second line represents 
the behavior of dMg{t)/dt in the instantaneous recycling approximation, where the returned 
(ejected) mass fraction, R^j, is approximately 0.7. 

The evolution of the gas mass is only weakly coupled to its chemical evolution. The 
returned mass, Mej{t), is almost independent of metallicity ( iKarakas fc Lattanzidl2007l ) and 
the stellar MS lifetimes are only weakly dependent on metallicity for values above 0.001. So 
with reasonable accuracy, the evolution of the gas mass can (but does not need to) be solved 
independently of the evolution of the dust and metals in the ISM. 



2.4. Integral Solutions 



When the evolution of the ISM gas is decoupled from its chemical evolution, eq. ([6]) can 
be readily solved with the aid of an integration factor, to yield a convenient functional form 
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for the evolution of the dust or any element in the ISM. 
Changing variables (jAudouze fc Tinsley 1976 ): 



dt 



dt 



dt 



eq. 



can be rewritten as: 



dZdjt) 
dt 



+ Zd{t) Fit) = Git) 



(16) 



(17) 



where Fit) and Git) are dimensional functions (units of time ^), given by: 



Fit) 



and 



m 

M,it) 
Git) 



1 + 



+ 



d,mw 



V^(t) dt 



dMgit) ^ 1 



dMgit) 

dr~ 



out- 



Mgit) 



Sit) + ZT^t) 



dMgit) 
dt 



inf 



(19) 



Note that the outflow term is included in Fit) because it is proportional to Zdit), whereas 
the infall term is included in the expression for Git), since it is independent of the evolving 
metallicity of the galaxy. Equation ( 1T7|) can be solved with the aid of an integration factor 
A(t), which is defined as: 



A(t) = exp 



Fit)dt 



givmg, 



Zait) = A(t)-i / A(t') Git') dt' + C 
Jo 



(20) 



(21) 



where C is an integration constant. 

Given Mgit), the mass of dust is simply given by: 



Mdit) = Zait) Mgit) 



(22) 



Finally, equation ([6]) can also be used to describe the evolution of the metallicity, Za, 
of any stable or unstable element A by the formal substitution of with M^, and by using 
the appropriate yields in eq. (fTTl) . For stable elements the parameter (md)^^ in eq. (fTSl) 
should be set to zero, and for an unstable element, the dust lifetime, r^, in eq. ([6]) must be 
substituted by its radioactive decay time. 

The solutions, as written here, are easily calculated for yields at a fixed metallicity, 
when the source function Sit) has a simple time dependency. Taking the changes of Sit) 
with metallicity into account can be achieved in a straightforward way. Yields are usually 
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calculated for a finite grid of metallicities {Zo{j),j = 1, A^}. One can therefore calculate the 
times {tj,j = 1,N} when the solutions derived with stellar yields at the fixed metallicities, 
M[Zo(j),t], reach the corresponding metallicities {Zo(j),j = 1,A^}. The final evolution of 
the mass, M{t), of the elements or the dust can then be written as a sum: 

N 

M{t) = J2 w,{t)M[ZoU),t] (23) 
1 

where {wj{t),j = 1,N} are weight functions chosen to ensure the continuity of M{t) across 
the grid of times {tj}. 

3. STELLAR SOURCES OF INTERSTELLAR DUST 

In the following we examine the yields and relative importance of the dust formation 
in SNe, WR, and AGB stars. In all our calculations we adopt the previously described 
mass- heavy IMF. 

3.1. Supernovae 

Consider a constant SN rate, Rsn- The dust production rate is simply given by 5*0 = 
Yd,sn X Rsn- Without grain destruction and with a SN rate of 20 yr~^, SN need to condense 
only 0.03 Mq of dust over a period of ~ 500 Myr in order to produce an observed dust mass 
of ~ 3 X 10^ Mq. When grain destruction by SN is taken into account, the amount of dust 
produced by SNe is given by: Ma ~ Yd^snRsN^d, giving the trivial solution: Y^^sn ~ ("^d), 
that is, in a steady state SNe must produce a dust mass that is equal to the amount they 
destroy during their lifetime. However, the required value of Y^^sn can be lower before the 
system reaches a steady state. For example, the first SNe expand in a dust-free medium and 
therefore are net producers of interstellar dust. Detailed chemical evolution and population 
synthesis models are therefore needed to fit the observational constraints imposed by the 
inferred stellar, gas, and dust masses, as well as the stellar and radiative output of the 
galaxy. 

Simple analytical models have shown that even in the absence of grain destruction, SNe 
must produce at least ~ 0.3 Mq of dust to account for the observed dust-to-gas mass ratio 
in J1148-I-5251 (Dwek et al. 2007; Figure 8). The reason that this number is not lower and 
equal to the value of ~ 0.03 Mq derived in our simple estimate above, stems from the fact 
that in these models the SFR is proportional to ^i'^sj where T^gas is the mass surface density 
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of the gas (iKennicuttl Il998bl ) . T^gas evolves with time, so that the effective timespan over 



which star-formation takes place is therefore much shorter than the age of the galaxy. It 
is determined by either the buildup of the ISM in infall models, or by the depletion of the 
ISM in closed box models. For Jl 148+5251 the situation is exacerbated by the fact that the 
IR-inferred SFR is ~ 3000 yr~^. This SFR can therefore only be sustained for about 
20 Myr, if the reservoir of gas is comparable to the dynamical mass or the CO mass of the 
galaxy. 

Do SNe actually produce this much dust? The most detailed information on the amount 
of dust formed in SNe comes from IR observations of the young supernova remnant of 
Cas A. Spitzer observat ions of the remnant have reveal ed the presence of ^ . 02 — 0.04 Mp, 
of ho t jRho et al.l 120081 ) and ~ 0.08 M© of cool dust dSibthorpe eTaPboOoi ISarlow et al" 



2010l ). Furthermore, detailed calculations, using a chemical kinetic approach to follow the 
transformation of gas phase molecules into small clusters of dust precursors, show that only 
~ 0.1 — 0.15 Mf^ of dust is created in the explosive ejecta of a Population III 20 Mq SN 
( ICherchneff &: Dwekll2010l ). Such yield may be typical of SNe of similar masses, regardless of 
the initial metallicity of the progenitor star. This yield is still not sufficient to account for the 
mass of dust observed in Jl 148+5251 except when grain destruction is almost nonexistent 
in that galaxy. Such special scenario is discussed in Section 4.1. 



3.2. Wolf-Rayet Stars 



Of the stars that become CCSNe, only those with masses above m^r ~ 40 Mq will 
become WR stars, a stage during they will experience extensive mass loss. The formation 
of carbon dust becomes feasible during the latest phases of this mass loss phase, when the 
stellar surf ace becomes car bon enriched, and the star evolves into a carbon-rich WR star 
(WC star) (ICrowtherl 119971 ). For an IMF with a Salpeter slope of a = 2.35, the fraction of 
CCSNe that become WR stars is: 



fu 



l-a 



0.11 



(24) 



So WR stars have to produce about 10 times more dust than SNe to account for the dust 
mass in J1148+5251. There is observational evidence that most of WR stars are part of a 
binary system with an OB companion and that the dust formation locus is located in the 
region where the two winds collide as exemplified by the archetypical systems WR104 and 
WR140. For the Pinwh e el Ne bula WR 104, a dust mass loss of ~ 8 x 10^^ Mq yr"^ is 



derived by iHarries et al.1 (120041). A typical 60 M (7^WR star is characterized by a WC phase 
lasting about 2 x 10^ yr (jPrantzos et al.lll986l ). We can thus approximate the amount of 
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dust formed by WR stars as ~ 0.2 Mq over their lifetime. This number is of the same order 
of magnitude and not 10 times larger than the dust mass formed in the ejecta of a CCSN. 
Therefore, WR stars are certainly contributing to the dust budget of Jl 148+5251 but no 
more than 10% of all dust makers. 



3.3. AGB Stars 



In the M ilky Way, a significant fraction of dust is produced in AGB stars f lDweklll998 
Tielenslll998l ). Dependent of their evolutionary stage on the AGB, red giants have either a 
oxygen- rich stellar photosphere characterized by a C/0 ratio < 1, or a carbon-rich photo- 
sphere resulting from thermal pulses and 3rd dredge- up and characterized by a C/0 ratio 
> 1. The former stars form silicate and metal oxide dust whereas the latter essentially form 
carbon dust and silicon carbides. At present, there exist no satisfactory models explaining 
the chemical formation of dust in AGBs and the impact of dust on accelerating the outflow 
through radiative pressure. Dynamical mode ls often use the classica l nucleation theory to 
describe dust synthesis and derive dust yields ( iFerrarotti fc Gailll2006l ). an approach that has 



been proved to be inadeq uate when applied to circumstellar outflows flDonn fc Nuthlll985 
Cherchneff fc Dweklboiol ). 



For the purpose of the present study, we ignore the chemical and physical complexity of 
the dust formation processe s and a ssum e a condensation efficiency of unity to calculate the 
dust yields as described by iDwekl (119981 ). The calculations assume that carbon dust forms 
when the C/0 number ratio exceeds 1, and that silicate dust forms when this ratio is less 
than unity. This approach provides upper limits for the dust mass. For illustrative purposes. 
Figure [T] (left panel) depicts the mass of dust produced during the AGB phase from 1 -8 M p, 
stars with an initial metallicity of 0.008. Yields were taken from lKarakas fc Lattanzid (120071 ). 
Also shown in the figure are the carbon and oxygen yields (in units of Mq and divided by 
their atomic mass), to illustrate the range of stellar masses that produce ei ther silicate or 



carbon dust. Yields similar to t hose presented in t his fig ure were adopted by lValiante et al. 



(120091 ) based on the models of IFerrarotti &: Gaill (120061 ). These models give different dust 
compositions, a distinction which is not important for the purpose of this paper. 

The right panel of the figure shows the IMF-weighted yield of dust from these stars, 
where the IMF is characterized by a power-law in mass with a Salpeter index of 2.35. The 
figure shows that he most efficient dust producers are ~ 3 Mq stars, which have a MS lifetime 
of about 350 Myr. For AGB stars to be significant contributors to the reservoir of dust in a 
galaxy, its age has to be significantly older than ~ 400 yr. 



10-' r 



10-' r 
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Fig. 1. — Left panel: The dust yield from AGB stars with an initial metallicity of 0.008 calculated using 



the yields of iKarakas fc Lattanzid (|2007( ). The red and blue lines depict the carbon and oxygen yields from 
these stars, divided by their atomic number. Carbon dust (hatched orange regions) is formed when the C/0 
ratio > 1, whereas silicate dust (hatched green region) is formed when that ratio is < 1. Right panel: The 
IMF-weighted yield of dust from these stars. The ^ 3 Mq stars, which have a MS lifetime of ~ 350 Myr 
(see Tabic 1) arc the dominant contributors to the production of dust. 



3.4. Relative Importance 



The rate of dust production by SNe, WR stars, and AGB stars is given by: ip x Yd^sn/iTT'i<, 
0.11 X ip X Yd^sn/iTT'-k, and ip x Y^^agb/ {^) , respectively. Taking 0.15, 0.2, and 0.04 Mq as the 
respective yields of dust in these sources, we get that in a steady-state, their relative dust 
production rate is: 1.0 : 0.11 : 4.6. AGB stars are therefore the most important sources 
of interstellar dust. However, this steady state is not immediately realized because of the 
delayed injection of AGB-condensed dust into the ISM. In the local group of galaxies, this 
delayed injection is manifested in the observed trend of the abundance of polycyclic aromatic 
hydrocarbons (PAHs) with galac tic metallicity, which is taken as a proxy for galactic age 
( lDwekll2005l : iGalliano et al.ll2008l ). The trend shows the existence of a threshold metallicity 
below which the abundance of PAH is very low. PAHs condense in the atmospheres of AGB 
stars, and therefore represent AGB-condensed dust. So at early times, the mass of dust 
i n galaxies is domina ted by SN-condensed dust. The dust evolution models presented in 
( iGalliano et al.l |2008| ) suggest that the steady state when AGB sources are the dominant 
dust producers may only be reached after about ~ (1 — 2) Gyr of stellar evolution. PAHs 
can also be produced by shock processing of AGB-condensed carbon dust. However, both 
mechanisms, stellar condensation and shocks, need AGB stars as the source of dust, so the 
PAH-metallicity trend manifests the delayed injection of dust by AGB stars, even if most of 
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the PAHs are formed in interstellar shocks. 



4. THE QUASAR J1148+5251 



4.1. The Evolution the ISM Gas and Dust Masses 



In the following we will apply our new mathematical formalism to calculate the evolution 
of the gas and dust in J1148+5251. For AGB stars to be an important source of interstellar 
dust, star formation must have commenc ed at least 400 Myr be fore the epoch of observations 
which is about 900 Myr. In their paper, IValiante et al.l ( 120091 ) adopted the SFH derived by 
Li07 from simulations of hierarchical galaxy mergers taking place at redshifts z > 6.5. The 
merger history of Li07 is not unique, and just one possible scenario leading to the formation 
of a supermassive BH at z ^ 6. We use it for illustrative purpose only, since it fulfills the 
necessary requirement for AGB stars to be significant sources of dust at high redshift. Their 
SFH, taken from their Figure 7 in Li07, is reproduced as a bold solid line in the left panel 
of Figure [51 

An important ingredient of chemical evolution models is the evolution of the gas mass. 
We adopted two different prescription for calculating this quantity. In the first, we adopted 
a closed box model (no infall or outflow) with an initial gas mass of 8 x 10^^ M©, and used 
eq. ( ITSll and the Li07 SFR to calculate its evolution. A value of R^j = 0.48 was chosen so 
that the stellar mass derived from the chemical evolution model would match that derived 
from the population synthesis model PEGASE (see Section 4.2). In the second approach we 
assumed that the SFR derived from the galaxy simulation follows the Schmidt-Kennicutt 
law in which ip{t) = /IM^^^(t), where A is a proportionality constant and a measure of the 
star formation efficiency in the galaxy. The value of A was chosen to reproduce the model's 
SFR of 60 Mq yr-i at t = 900 Myr with a gas mass of ~ 5 x 10^° Mq. 

In the calc ulations we adopted an IM F-averaged dust yield of ld,sn = 0.15 Mq, the 
yield derived by ICherchneff fc Dwekl f l2010l ) for a 20 Men p r imord ial core-collapse supernova. 
The yiel d of dust in AGB stars was c alculated as in iDwekl (jl998[ ) using the stellar elemental 
yields of iKarakas fc Lattanzid ( 120071 ). The total dust yield from AGB stars is only weakly 
dependent on metallicity, so for simplicity we adopted AGB yields at a constant metallicity 



oi Z = 0.008. The stellar IMF was chosen to be a power law: 
100 Mq, with a spectral index of a = 2.35. 



between 1 and 



Figure [H depicts the evolution of dust and gas masses for the two adopted evolutionary 
models for Mg{t): the closed box model (left panel), and the model described by the Schmidt- 
Kennicutt SF law (right panel). The vertical line at t = 900 Myr is the observed dust 
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redshift redshift 
18 14 12 10 8 7 6.5 18 14 12 10 8 7 6.5 




200 400 600 800 1000 200 400 600 800 1000 

time (Myr) time (Myr) 

Fig. 2. — The evolution of the mass of the SN-condensed dust, AGB dust, total dust mass, and the mass 
of the ISM gas in J1148+5251, for the SF scenario of Li07, depicted here in Fig. ([5]). Left panel: Results 
for the closed box model. Right panel: Results for the Schmidt-Kennicutt SF law. The vertical line at 
t — 900 Myr is the observationally-infcrrcd dust mass. 



abundance, t he length of the b ar reflecting the range of dust masses for the different dust 



composit ions (iDwek et a" 



20071) . The results of our calculations are in good agreement with 



those of IValiante et al.l ( l2009l ). Not surprising, the detailed evolution differs from theirs, 



since it depends on their adopted evolutionary history of the gas mass, the grain destruction 
efficiency, and AGB yields, all quantities that differ somewhat from the ones used in this 
work. However, the main results, the final masses produced by SN- and AGB-condensed 
dust, are very similar. The figures show that AGB stars can readily produce the observed 
amount of dust with the proposed SFR, and that SNe contribute less than 20% of the total 
dust mass at t = 900 Myr. They can therefore only be important dust sources if their 
IMF-averag e d yie ld is ~ 1 Mq, confirming the results of the simpler analytical model of 



DweketaL 



(120071 ). 



Figure [2] also shows that for the Schmidt-Kennicutt model there exist two epochs, oc- 
curring ~ 100 and ~ 400 Myr after the commencement of SF at z ^ 15, during which SNe 
produce over ~ 10*^ Mq of dust. It therefore seems that in a SF scenario consisting of either 
an extremely young or an intermediate-age burst of stars, SNe can produce the observed 
amount of dust, even with an average dust yield of ~ 0.15 Mq, without resorting to the need 
of AGB stars. 

At the time of the intermediate burst, the SFR is over 10^ Mq yr~^. Compared to the 
younger burst, this burst requires a significantly larger SFR to produce a similar amount 
of dust, a consequence of the effect of grain destruction. The corresponding bolometric 
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luminosity is over 2 x 10^^ Lq. Both, the SFR and the luminosity are in excess of the 
observational constraints, ruling out this SN scenario as a viable source of dust. 

The extremely young burst will produce a total bolometric luminosity of ~ 5 x 10^'^ Lq, 
similar to the observed luminosity of ~ 10^^ Lq, and a stellar mass of ~ 10^^ Mq (see 
Figure [3]). This is a viable scenario for the formation of dust in J1148+5251, although it 
requires the formation of supermassive BH on a very short timescale, and has implication 
on the comoving number and luminosity densities of such objects in the early universe (see 
Section 6.2 below). 



4.2. The Evolution the Stellar Luminosity and Mass 



Figure [3] depicts the evolution of the stellar masses and remnants (left panel) and the 
bolometric stellar luminosity (right panel) of J1148+5251. Calculations were done for the 
hierarchical SF history depicted in figure using the PEGASE population synthesis code 
(IFioc fc Rocca-Volmerangd 119971 ). Similar figures with identical results were presented by 
Li07 (Figures 11-12 in their paper). The merger history was proposed by the m to account 
for th e presence of a supermassive black hole (BH) of mass ~ 3 x 10^ Mq (jWillott et al. 
20031 ). Assuming that the relationship between the BH mass and its spheroidal component 
is s till valid at redshifts > 6, giv es a stellar bulge mass of ~ few x 10^^ Mq [e.g. review 
by (IKormendy fc Gebhardtl 120011 )]. While the stellar mass derived here agrees with that 
estimate, it is in excess of the dynamica l mass of ~ (5.0 ± 2. 5) x 10^° M(7) that is enclosed 



within a 2.5 kpc radius of J1148+5251 (IWalter et al.ll2004l ). IWalter et all (120041 ) proposed 
several solutions to this discrepancy, including an overestimate in the mass of the BH, and 
a breakdown in the assumption that the CO gas is gravitationally bound. Also, there may 
be evidence that at high redshifts BH s may grow very r apidly without the corresponding 
increase in the mass of the host galaxy (jShields et al.ll2006l ). It is therefore premature to rule 
out the AGB dust formation scenario while the origin of this discrepancy is still unresolved. 



5. THE DEPENDENCE OF DUST MASS ON THE 
STAR FORMATION HISTORY 

As pointed out earlier, the Li07 SF scenario for the formation of J 1148+5251 is not 
unique. We therefore explore the dependence of the total dust mass on the SFH of the 
galaxy. We first consider the production and survival of dust by single short bursts of star 
formation. Figure H] depicts the evolution of the mass of dust produced by AGB stars that 
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redshift redshift 
18 14 12 10 8 7 6.5 18 14 12 10 8 7 6.5 
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time (Myr) time (Myr) 

Fig. 3. — Left panel: The evolution of the stehar mass content (left panel) and luminosity (right panel) 
calculated with the population synthesis code PEGASE using the hierarchical SFH depicted in Figure [5l 



were born at time t = in a 100 Myr burst of star formation characterized by a constant 
SFR of 1000 Mq yr~^. The bold sohd hne depicts the net production of dust by AGB stars 
when no grain destruction occurs. It takes about ~ 50 Myr for the first AGB stars to evolve 
off the main sequence. The thinner lines depict the dust evolution in the presence of grain 
destruction. Grain destruction was assumed to be constant and to commence only after 
the cessation of the burst. After 600 Myr the burst releases ~ 1 x 10^ Mq of dust if no 
destruction has taken place in the intervening period since the end of the burst and the time 
of observation. The surviving mass of dust is significantly smaller when grain destruction 
is taken into account. A young burst produces after 200 Myr the least amount of dust, 
about 10^ Mq with no grain destruction, but is also least affected by grain destruction. The 
figure shows that a single burst of star formation (with no subsequent grain destruction), 
for example burst 1 with an intensity of ~ 3000 Mq yr~^, or burst 4 with an intensity of 
10,000 Myr, can produce the observed ~ 3 x 10^ Mq of dust in J1148+5251. 

Deconstructing the complex SFH of a galaxy undergoing a series of hierarchical mergers 
can be very illustrative for understanding the origin of dust in such system s . In F igure [5] 



(left panel) we approximate the SFH of J1148+5251 used by IValiante et all (120091 ) by five 
discrete bursts of star formation. The right panel depicts the cumulative contribution of 
each of the bursts to the total mass of dust in the galaxy. The oldest burst (number 1) 
releases the most amount of dust since all stars above 2.1 Mq contributed to its production. 
However, grain destruction by the subsequent bursts reduces the mass of dust that survives 
until 900 Myr to only ~ 2 x 10^ Mq. Almost all of the surviving dust was produced by the 
third burst. It has a main sequence turnoff mass of ~ 3 Mq, but also the largest SFR, and 
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Fig. 4. — The evolution of the mass of dust released by AGB stars in a burst of star formation that 
commenced at < = as a function of time. The burst duration is 100 Myr, and its intensity is 1000 Mq yr^^. 
The top bold line represents the evolution of the dust mass without any grain destruction. Other lines depict 
the evolution of the dust when it is destroyed by a constant ongoing destruction mechanism, assumed to 
have started after the end of the burst. The lines are marked by the dust lifetime in units of Myr. The 
vertical dashed lines mark the age of the burst at epochs that correspond to the burst ages in Fig. [5] 
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it is followed by a period of relatively low star formation activity and grain destruction. 

The formation of about ~ 3 x 10*^ Mq of dust at 2; > 6 can therefore be achieved 
for different SFHs. As illustrated in Figure HI a single burst of star formation with the 
appropriate intensity can produce the required amount of dust, since no grain destruction 
takes place after the cessation of all star formation activities in the galaxy. Massive amounts 
of dust can also be produced by complex SFHs such as that presented in Figure |5l provided 
that the epochs of intense starburst acivity are followed by a period with a very low rate, 
< 70Mq yr~^, of star formation, so that the dust lifetime will be over 100 Myr. The 
discriminating factor between the different star formation histories capable of producing 
large amount of dust are the cumulative products, such as the metals or stellar masses, 
resulting from the star formation activities. 
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Fig. 5. — Left panel: The star formation rate of a galaxy going through hierarchical merger episodes 



(Li07, see also Figure 1 in (|Li et al.ll2008l )) shown as a bold red line, can be approximated by a series of 5 
discrete 100 Myr duration bursts (shaded bars). Right panel: The contribution of each of the five bursts 
depicted in the left panel to the total dust mass observed at the galactic age of 900 Myr. Most of the dust 
is that created by the third burst which had the largest star formation rate, and is followed by a period 
of relatively low rates of star formation and grain destruction. The figure illustrates the sensitivity of the 
surviving dust mass to the SFH of the galaxy. 
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6. OBSERVATIONAL DISCRIMINATION BETWEEN 
THE SN AND AGB SCENARIOS 

We have shown that not only AGB stars, but also SNe can produce the inferred dust 
mass in Jl 148+5251, even with observed SN dust yields inferred from observations of Cas A. 
The SN scenario requires a special SF event, characterized by an intense burst of star for- 
mation that commenced just prior to the observations. The short duration of the burst can 
insure that grain destruction did not have any time to significantly erode the newly-formed 
dust. This is clearly an idealized scenario, since s ome grain destructio n will take place when 



the reverse shock travels through the SN ejecta ( INozawa et al.l 120101 ). In the following we 



discuss several observational effects that can discriminate between the two scenarios. 



6.1. The Spectral Energy Distribution 



An important distinction between the SN and AGB scenarios for the origin of dust in 
J1148+5251 is the origin of its SED. In the SN scenario, most of the UV-optical and all the 
far-IR emission originate from starlight and starlight-heated dust. In the AGB scenario the 
situation is more complicated, and the origin of the SED depends on the fractions of the 
AGN and starburst luminosities that are absorbed and re-emitted by the dust. The total 
calculated stellar luminosity a.t z = 6.4 {t ~ 900 Myr) is 1.6 x 10^^ Lq (see Figure[3]and Table 
|2]), about equal to the total far-infrared (FIR) luminosity, which is ~ (2.0 ± 0.5) x 10^^ Lq. 
So, in principle, all the FIR emission could arise from dust heated by the stellar radiation. 
In this extreme case, all the observed UV to near-IR luminosity should originate from the 
AGN. Figure [6] (left panel) depicts the stellar SED calculated for the hierarchical merger 
SFR at the galactic age of 900 Myr. In the absence of any discriminating line emission, 
can this extreme case be ruled out on the basis of the 0.1 to 1 /im continuum emission? 
This question is addressed in the right panel of the figure which depicts the probability 
distribu tion of photornetric spectral indices ap derived from the observed magnitudes of 



quasars 



Richards et al.l (120031 ). The spectral index of the 0.1-1.0 /xm spectrum of J1148+5251 



is ~ —0.35 (left panel) which corresponds to a value of Op ~ —1.65. This value is near the 
median value of -1.6. The optical to near-IR spectrum of J1148+5251 has therefore an equally 
high probability of being entirely that of an AGN as that of a starburst galaxy. The nature 
of the SED can therefore only be determined with the detection of mid-IR fine structure 
lines from highly ionized spe cies such as [Ne V]14.3 and 24.3 jum w hich are unambiguous 
diagnostics of AGN activity flSatvapal et allbood : budik et al.ll2009h . 
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Fig. 6. — Left panel: Comparison of the stellar SED calculated for the hierarchical merger SFR at 
2;=6.4 (shaded area) with the observations. In the AGB most of the optical to near-IR flux is emitted by 
the AGN. In contrast, in the SN scenario for the origin of the dust, most of this emission originates from 



stars [see Fig 11 in (jPwek et al.ll2007l )]. Right panel: The probability distribution of photometric spectral 



indices ap derived from the observed magnitudes of quasars [Richards et al.l (|2003l) . The spectral index of 
a,y w —0.35 (left panel) corresponds to a value of ap ?s —1.65. The figure illustrates the fact that there is a 
high probability that the observed 0.1-1 /umflux could be entirely attributed to the emission from an AGN. 



6.2. The Rarity of Dusty Hyperlumimous IR Galaxies at High Redshifts 

The AGB scenario for the forma tion of dust in J1148+5251 requires the early onset of 
star formation. IValiante et al.l (120091) used the SFR derived in the Li07 hierarchical mergers 
model as a plausible description of the SFH in J1148+5251. Figure E] shows that this scenario 
requires the assembly of galaxies with stellar masses of M^, ^ 7 x 10^^ Mq by z ~ 9. 



The associated mass of the dark matter (DM) halo is: Mh ~ (Qm/^b) {M^ + Mg) ^ 4 x 
10^^ Mq. The com oving number density of such halos derived fro m the Press- Schechter 
(PS) approximation (IPress fc Schechterlll974h is about 10"^^ Mpc"^ (Istiavellilbood : Ihi et al. 
20071 ). This number is sig nificantly smalle r that the ~ 10^^ Mpc~'^ comoving volume density 
of the 2 ~ 6 SPS S QSOs (|Fan et al.ll2003r) which have optical to near-IR luminosities similar 



to J1148-I-5251 (j Jiang et al.ll2010l ). The apparent three orders of magnitude discrepancy 
between the observed number density of J1148+5251-like QSOs and that calculated using 
the PS approximation can probably be attributed to one or more of the following: (1) the 
breakdown in the accuracy of the formalism at high reds hifts and large hal o masses (Li07, 
and references therein); (2) the observations presented by I Jiang et al.l (120101 ) only cover the 
optical to near-IR wavelengths, which only includes a hot dust component. Their QSO 
sample may therefore not contain a massive amount of colder dust, and may therefore not 
be representative of the population of J1148+5251-type objects at far-IR wavelengths; and 
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(3) the Li07 merger scenario is not the only scenario that can lead to the formation of the 
supermassive black hole, even if the BH-halo mass correlation holds at these redshifts. 

In the SN scenario, if we assume that the dust formed in an extremely young burst, the 
mass of the dark halo associated with an ^ 10^^ Mq galaxy is about 6 x 10^^ Mq. The 
comoving number density of such halos is about 10~^ Mpc~^, consistent with the number 
density of galaxies with simila r stellar masses d etected in the Great Observatories Origins 
Deep Survey (GOODS) fields (IStark et al.ll2009l ). If all these galaxies had luminosities sim- 
ilar to J1148-I-5251, then the comoving bolometric luminosity density at 2; ~ 6 should be 
about 10^ Lq Mpc~^. This lumin osity density is soni e what larger than the total luminosity 
density presented in Figure 18 of iFranceschini et al.l ( 120101 ). when extrapolated to z ~ 6. 
However, galaxy count models suggest that galaxies with luminosities above ~ 2 x 10^^ Lq 
constitute only a smal l a fraction, about 10~'^, o f the total number of galaxies at high redshift 
(jChary &: Elbag|200ll : IFranceschini et al.ll2010l ). This gi ves a comoving n umber density of 
J1148+5251-type objects of ~ 10~^ Mpc~^, similar to the iFan et al.l ( l2003l ) observations. So 
in the SN scenario a significant fraction of the 2; ~ 6 quasars should have massive amounts 
of dust. This result is not yet supported by current observations, which show that a signif- 
icant fraction of these objects are dust-free. However, among the ^ 10^^ M© galaxies, 
J1148-I-5251 is an extremely rare object, consistent with the inferred low extinction and 
metallicities in z ~ 7 galaxies detected i n the CDF-South GO ODS field with the WFC3/IR 
camera on the Hubble Space Telescope (IBouwens et al.ll2010l ). The rarity of Jl 148+5251 in 
the SN scenario should not be surprising considering the very contrived SFH required to 
form the inferred mass of dust in this object. 



6.3. The Star Formation Efficiency 



Observations of star forming galaxies revealed an empirical relation between a galaxy's 
global SFR and i ts total gas mass. For a Salpter IMF, this relation can be written as 
KennicuttI Jl998bh : 



^SFR 



{Mq yr ^ kpc 



(2.5 ±0.7) X 10~^ Sg(M0pc 



-2\1.4±0.15 



(25) 



VLA observations of J1148+5 251 show that the C O gas is distributed over an extended 
region with an area of ~ 4 kpc^ (IWalter et al.ll2004j ). Adopting a total gas mass of 5 x 
10^° Mq, and assuming that it occupies the same area as the CO gas gives an expected SFR 
of ~ 560 ± 240 Mq yr~^. For the mass heavy IMF adopted here, the corresponding SFR is 
~ 220 ± 90 Mq yr^^. This value is larger than the SFR at the epoch of observations in the 
AGB scenario, suggesting that the SF efficiency at that time is about ~ 3 times lower than 
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that in normal star forming galaxies, which varies between 2 and 10% f lKennicuttl Il998bl ) . 
However, this instantaneous SF efficiency does not reflect its past value which, averaged over 
a fixed time interval, should scale as So the average SF efficiency in J1148+5251 should 
be similar to that in normal SF galaxies. 

Studies of the relation between the mass of dark matter (DM) halos and the luminosity of 
their host galaxies show an evolution ary trend between the galaxies' luminosity function and 
their halo masses (ITrenti et al.ll2010l ). This trend can be translated into a relation between 
their SF efficiency, r], and their halo mass. For masses above 10^^ Mq, the SF efficiency 



depen ds only mildly on halo mass, and is given by ?7 ~ 0.06 (M/i/10 Mq) fiTrenti et al. 



2010l ). For the inferred halo mass of 4 x 10^^ Mq in the AGB scenario, r] f» 0.18. So the 
predicted SF efficiency is somewhat higher than the inferred average in the AGB scenario, 
which can however be attributed to expected large variations in individual objects. 

In the SN scenario, all the far-IR luminosity is powered by star formation, suggesting 
a SFR in excess of ~ 3000 Mq yr~^. This implies a very high SF efficiency for the given 
mass of the galaxy. At this rate the galaxy will deplete its reservoir of gas in about 10 Myr, 
suggesting that we are just witnessing the onset of a large burst of SF in this object, which 
will be rapidly quenched by stellar feedback. 



7. ALTERNATIVE DUST SOURCES 

The still unresolved discrepancy between the predicted and inferred stellar masses in the 
AGB scenario, and the specialized SF scenario required to form the inferred Jl 148+5251 dust 
mass by SNe, raises the need to explore alternative, non-stellar, sources that may produce 
the large amount of dust in this object. In the following we consider the growth of dust in 
molecular clouds and their formation around AGN as such sources. 



7.1. Molecular Clouds 



Numerous observations and theoretical considerations support the notion that dust 
grains are processed and grow in dense molecular clouds (MCs). Observationally, extinc- 
tion measurements show that the value of Ry = Ay/EjB — V) , an in dicator of grain 



size, increases with hydrogen column density [e.g. iKandori et al.l (j2003[ ). and references 
therein]. Simila rly, interstellar depletion s of refractory elements increase with average gas 
column density f lSavage fc Sembachlll996t |Jenkinsll2009l ) and infrared spectroscopy of molec- 
ular clouds show the presence of ices, expected to have condensed out on interstellar grain 
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cores (Ivan Dishoeckl 1200 



)■ Theoreti cally, the discrepancy between the relatively short life- 



time of interstellar dust ( lJonesll2004f ) and the longer lifetime for their r eplenishment by SNe 



or AGB sta r s has been at t ributed to grain growth in molecular clouds (IDwek fc Scald Il980 



Dwek 1998; Tielens 1998; Zhukovska et al. 



20081 ). In particular, the possibility that rapid 



growth of dust in high-r edshift galaxies took place in molecular clouds was considered by 
Michalowski et al.l (120101 ) . Here we provide a more quantitative discussion of such possibility. 



Several conditions must be met for grain growth in clouds to be significant: (1) the 
presence of preexisting refractory grain cores onto which the accretion can take place; (2) 
the accretion time onto these cores must be shorter than the lifetime of the dense clouds in 
which this process is taking place. 

Let 

Pgr — "^gr ^ '^^gv 

be the mass density of grains in the MC, where their number 

density and rrigr their mass. Assuming that all grains have the same radius a, the growth 
rate of pgr due to accretion of a heavy element A with a molecular weight fiA and number 
density ua is given by: 



dp 



gr 



dt 



a Tia Pa ^gr v 



where a is the sticking coefficient in the collision, and v 
velocity of A at temperature T. 



(26) 

{SkT/npAY^^ is the mean thermal 



The accretion time, Tacc, is defined as follows: 



1 dPgr 
Pgr 



dt 



a 



a 



vra 



gr 



Apa 



Pa nAV 
Pa 



(27) 



riA 



Taking the mass of an oxygen atom to represent that of a colliding species, so that pa/^h 
16, and nAjnu^, ~ 5 x 10~^ (assuming solar abundances), we get that: 

18 / ^ \ . . _i /T-\ -1/2 



2 X 10 



Ta 



sec 



P 



a \ f 



(28) 



a \pm J \cm 

For a grain density of 3 g cm~^, radius of 0.1 /im, a cloud density of 10^ cm~^, and cloud 
temperature of 20 K, the accretion time becomes Tacc ~ 3 x 10^ yr. 

The lifetime, Tc, of the molecular clouds is given by: 



2 X 10 



10 



7 X 10^ yr 



(29) 



^ 3000 

where is the MC mass, and where we used the inferred H2 mass and SFR to calculate 
the average lifetime of the molecular gas. We used here the large SFR inferred from the 



- 24 - 



SN scenario to minimize the lifetime of the MC. Had we used the SFR from Figure [5] the 
MC hfetime would have been higher by a factor of ~ 50. Within the uncertainties of the 
chosen values of the various parameters, the cloud lifetime is not significantly shorter than 
the accretion time. Consequently grain growth in the dense ISM can play an important role 
in determining the observed dust mass. 



The accreted material may be in the form of organic refractory material ( iGreenberg et al. 



19951 ). The main observational constraints, extinction, diffuse IR emission, interstellar abun- 
dances, i n the Milky Way can a l so be met with com posite dust models that incorporate such 
material Jli fc Greenberg|ll997l : IZubko et all2004h . If so, then the dust of J1148+5251 could 
also consist of such composite grains that could have a significant abundance in the diffuse 
ISM. In the Milky Way, most (60-75%) of t he IR emission is rad iated by dust in the diffuse 
ISM, with only 15-30% radiated from MCs ( ISodroski et al.lll997l ). Likewise, the IR emission 
from Jl 148+5251 could arise either from the dust in the MCs, or from dust in the diffuse 
ISM, including H II regions. The origin of the IR emission in this case can therefore only 
be determined by considering dust evolution models that incorporate a two phase medium, 
and take the physical processes and the cycling between the phases into account. 



7.2. "Smoking Quasars" 



The concept of "smoking quasars", that is, t he possibility that quasars may be pro- 
ducers of interstellar dust was first put forward by lElvis et al.l ( 12002| ). In his scenario, the 
broad emission line clo uds (BELCs), confined by the centrifugally driven outflowing wind 
( iKartje fc Konigllll996l ). go through the same density and temperature conditions that char- 
acterize the dusty winds in stellar outflows. It is therefore naturally to postulate that if the 
conditions in the stellar winds are ripe to form d ust, a similar process will be taking place 
in the BELCs. This idea was further discussed by lMaiolino et al.l ( l2006l ). 



The main issue, not address ed in these papers is whether the BELCs were dusty to start 
with (lElitzur fc Shlosmarul2006[ ). If so, then the IR emission observed in the vicinity of the 
AGN is not the signature of radiation from newly-formed dust, but that of reprocessed dust. 
This dust could have formed previously in SNe or dusty stellar winds. If so, quasars cannot 
be considered as net producers of interstellar dust. The issue of quasars as dust sources can 
therefore only be settled once we have a clearer understanding the origin of the BELCs, and 
of the dust formation process in these objects. 
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8. DISCUSSION AND SUMMARY 

Using our new integral solutions for the chemical evolution of the ISM gas and its metal 
and dust contents, we examined the origin of the large amount of dust discovered in the 
high— 2; quasar J1148+5251. We have confirmed that AGB stars can produce the required 
amount of dust in a merger scenario in which the assembled galaxy undergoes a period of 
intense star formation at 2; ~ 9, followed by a period of significantly lower stellar activity. In 
this AGB scenario, an average 20 Mq supernova needs to make at least ~ 1 Mq of dust in 
order to be a viable source of dust at this high redshift. However, SNe can yet produce the 
required amount of dust with a significantly lower dust yield, provided the galaxy underwent 
a very special SFH. In the following we summarize the main results of our paper, and discuss 
future observations and studies needed to distinguish between the two scenarios in order to 
determine the origin and nature of the dust in Jl 148+5251 and similar objects in the early 
universe. 

1. The AGB scenario - A galaxy must be at least 400 Myr old for AGB stars to have 
made a significant contribution to its reservoir of dust, a conditioned fulfilled by the 
merger scenario for the formation of J1148+5251. However, this SFH is not unique, 
and we outlined several alternative scenarios for the formation of large amounts of dust 
with AGB stars. At any given epoch, the intensity of the starburst must be at least 
2000 Mq yr~^ in order to produce the inferred amount of dust, and the star formation 
rate must have dropped to less than ~ 50 Mq yr~^ thereafter so that the dust that is 
injected at the epoch of observations will not be significantly destroyed. 

2. The SN scenario - In general, an average SN must condense about 1 Mq of dust 
in order to produce the observed amount of dust in J1148+5251. Our current state 
of knowledge, observationally and theoretically, suggest that the explosive ejecta of a 
progenitor ~ 20 Mq star produces only about 0.1 — 0.15 Mq of dust. Such progenitors 
are the most common metal producing SNe, and fall short of meeting the production 
requirement. However, our knowledge is based on limited number of observations and 
calculations, so it may change in the future. Similarly to the alternative AGB scenario, 
SNe can produce the observed amount of dust with a dust yield of ~ 0.15 Mq in a 
single short-duration and intense burst of star formation, followed by a period of very 
low star forming activity. Such scenarios must be very rare in order to explain the low 
comoving number density of hyperluminous IR objects at high redshifts. 

3. The stellar mass and gas content of J1148+5251 - The stellar mass offers an 
important integral constraint on the duration and intensity of the burst of star for- 
mation. These place integral constraints on the total amount of dust that could have 
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formed in the quasar. In the SN scenario ^ 10^^ Mq of stars are formed around 
2; ~ 7. So if dust formation was very efficient in these objects, the SN scenario would 
produce more dusty hyperluminous IR stars than currently observed at high redshifts. 
The formation of SN dust must therefore be a very rare event in these objects, limited 
by the dust formation and grain destruction efficiencies in these objects. The AGB 
scenario requires the formation of 7 x 10^^ Mq at redshifts of ~ 8 — 9, suggest- 

ing that J 1148+5251 is a very rare event, but perhaps consistent with the number of 
J1148+5251-like QSOs detected at 2; ~ 6. If so, a significant fraction of 2; ~ 6 quasars 
should have massive amount of dust, a prediction that is yet unsupported by current 
observations. Furthermore, the stellar mass produced in the AGB scenario is much 
larger that the inferred dynamical mass of the galaxy. Resolving this discrepancy is a 
crucial step in determining the origin of the dust in Jl 148+5251. 

4. The nature of the optical to near-IR SED of J 1148+5251 - In the SN scenario 
all the emission in this wavelength region is produced by stars, whereas in the AGB 
scenario all this emission is produced by the AGN. In the absence of a detailed spec- 
trum, we used the slope of the SED to discriminate between the two scenarios. The 
results show that the slope of the SED is typical of quasars, and also consistent with 
that produced by an extincted burst of star formation. The detection of mid-IR fine 
structure lines that serve as unambiguous diagnostics of AGN or starburst activity will 
be extremely useful for distinguishing between the two scenarios. 

5. The efficiency of grain destruction in the ISM - Our current state of knowledge is 
predominantly based on the simulations that were performed for a homogenous ISM. 
These calculations need to be expanded for a 2-phase ISM, taking also the spatial 
distribution of dust sources and sinks into account. For example, one can easily envision 
an intense starburst or spatially isolated SNe producing large amount of dust in one 
location in the galaxy, followed sometime later by another physically isolated starburst 
that will have little effect on the destruction of the dust created in the first. Spatially 
isolated starburst or SNe may therefore increase the currently estimated short lifetime 
of the dust in the ISM. 

6. Dust growth in the ISM - Accretion in molecular clouds must play an important role 
in explaining the origin of the Galactic pattern of interstellar depletions. A necessary 
condition for this process to be important is that the lifetime of molecular clouds 
be larger than the timescale for grain growth by accretion. This condition is met in 
J1148+5251, rendering it a possible, but not necessary, source for the growth of dust 
in this object. Detailed information of the correlation of the IR emission with the 
different gas phases in Jl 148+5251 may resolve this ambiguity. 



-27- 



7. Dust formation in around AGN - The effectiveness of AGNs as sources of newly- 
formed dust remains unclear. Tlie origin of the dense broad-line-emitting clouds, and 
whether they were initially dusty so that the AGN merely reformed pre-existing dust, 
are questions that need to be addressed to ascertain their role as dust sources in quasars. 
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Table 1. Main- sequence Lifetimes of AGB Progenitor Stars^ 



Mass 
(Me) 






Metallicity 






0.0001 


0.001 


0.008 


0.0100 


0.02 


1.0 


4998 


6774 


8634 


8846 


9516. 


1.5 


1474 


2022 


2415 


2446 


2528. 


2.0 


663 


914 


1047 


1054 


1062. 


2.25 


486 


670 


756 


759 


757. 


2.5 


371 


512 


570 


570 


565. 


2.75 


292 


404 


444 


443 


436. 


3.0 


236 


326 


355 


354 


346. 


3.25 


195 


270 


291 


290 


282. 


3.5 


164 


ZZ { 


z43 


Z41 


Z64:. 


3.75 


140 


194 


206 


204 


197. 


4.0 


121 


168 


177 


175 


168. 


4.5 


94 


129 


135 


134 


128. 


5.0 


75 


104 


107 


106 


100. 


5.5 


62 


85 


87 


86 


81. 


6.0 


52 


72 


73 


72 


68. 


6.5 


45 


61 


62 


61 


57. 


7.0 


39 


53 


53 


53 


49. 


7.5 


34 


47 


47 


46 


43. 


8.0 


31 


42 


41 


41 


38. 



^ Lifetimes in Myr, calcu l ated using the analytical 
formulae from Raiteri et al. ( 19961 ). 



Table 2. Energy Output from J1148+5251 



Wavelengths Luminosity {Lq) 

UV-near-IR (0.1-5 ^m) (5.7 ± 1.3) x lO^^ 

mid-IR (5-30 /im) (1.8 ± 0.2) x lO^^ 

far-IR (30-300 /im) (2.5 ± 0.5) x 10" 

bol (0.1-300 Aim) (10.0 ± 1.5) x lO^^ 



calculated stellar 1.6 x 10 



13 



^The total bolometric output at t = 900 Myr 
(see Figure [3]). 



